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Multi-Species entanglement, defined for a many-particle system as the entanglement between 
different species of particles, is shown to exist in the thermodynamic limit of the system size going 
to infinity. This macroscopic entanglement, as it can exhibit singular behavior, is capable of tracking 
quantum phase transitions. The entanglement between up and down spins has been analytically 
calculated for the one-dimensional Ising model in a transverse magnetic field. As the coupling 
strength is varied, the hrst derivative of the entanglement shows a jump discontinuity and the 
second derivative diverges near the quantum critical point. 
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Entanglement in a quantum state is a signature of 
quantum correlations between different parts of the sys- 
tem. Quantum entanglement, perceived as a resource for 
quantum communication and information processing, has 
emerged over the last few years as a major research area 
in various diverse fields such as physics, mathematics, 
chemistry, electrical engineering and computer science 
0, H|. The von Neumann entropy of a subsystem is a 
measure of the entanglement between the subsystem and 
the rest. Mostly, the entanglement between two spatially 
separate and distinct parts has been studied in various 
systems. In this article we will address the entanglement 
between two species of particles sharing the same physi- 
cal space. As we shall see below, that the thermodynamic 
limit exists for the two-species entanglement, i.e. the en- 
tropy per volume (size) tends to a finite constant, in the 
thermodynamic limit of the system size becoming macro- 
scopically large. Thus, the macroscopic two-species en- 
tanglement is a natural candidate for tracking quantum 
phase transitions, which will be explored in the context 
of the transverse-field Ising model of the spin systems, a 
prototypical system of quantum critical phenomena. 

Let us consider two species of hard-core particles A and 
B moving on a lattice of N sites, the corresponding num- 
bers of particles for the species being Na and Nb respec- 
tively. Local hilbert space for a site is four-dimensional, 
either A or B occupation, double occupation, and no oc- 
cupation. The basis states of the bipartite system can 
be chosen to be the direct products of the basis states 
of the individual parts: \u)a\v)b, where u (v) is a set of 
site locations occupied by A-type (B-type) particles. A 
general pure state of the combined AB-system is given 
by 



^2ip(u,v) \u) A 



\v)i 



(1) 



and the wave function amplitudes ip(u,v) determine all 
the properties of the system. The von Neumann en- 
tropy of the subsystem (of particles A), the measure 
of the entanglement between the two species is calcu- 
lated from the eigenvalues of the reduced density matrix 
Pa = Trs\'4>){'4>\- In general the computation of pa is 
quite involved, but it is easier for the case of exclusion 
and half filling. That is, the exclusion forbids double oc- 



cupancy of the two species at a given site, and the con- 
dition of half filling (the total number of particles equals 
the number of sites) implies that every site is occupied 
by either A or B particle. In this case, the exclusion and 
half filling conditions map all the physical states of this 
system to a spin system, each site occupied by a spin- 1/2. 
Since, each site is either occupied by one of the species, 
knowing the locations of A particles is enough to charac- 
terize the wave functions shown above. Here, v would be 
just the complement of u. Now, the state can be written 



U)A \U)b, 



(2) 



where for the state of B, the set u denotes site loca- 
tions not occupied by B particles. Now, the above bipar- 
tite state already has the desired Schmidt decomposition 
built in. The quantum entanglement between the two 
species is given by the von Neumann entropy, 



s A = -]T|v>WI 2 iog|v>WI 2 



(3) 



This is the desired simplification for the exclusion at 
half filling, with the wave functions themselves being the 
Schmidt numbers. The entanglement between the two 
species increases with the correlation present in the state. 
In the case of the two species being up and down spin 
particles of many-electron systems, a strongly-correlated 
state (with severe local constraint of no double occu- 
pancy) will have larger entanglement and uncorrelated 
spin states (for example metallic states) do not have any 
entanglement. In contrast, the usual entanglement mea- 
sures studied for distinct spatial blocks is maximized if 
the local constraints are relaxed 

The above von Neumann entropy is invariant if the ba- 
sis chosen for A or B particles is the momentum basis (a 
unitary transform of the position basis used here), but 
it is not invariant under a SU(2) unitary that can lead 
to hybrid particles. This should be contrasted with the 
usual entanglement in spin systems explored between two 
distinct spatial parts, which is invariant under a SU(2) at 
any site, but not invariant if the basis is changed to the 
momentum basis (as this will mix the two partitions). 
Since both species of particles access the full physical 
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space of states, independent of the basis, the entangle- 
ment between the two species is expected to be extensive 
with the system size. If the density of particles for both 
species is finite, there is a possibility of a macroscopic en- 
tanglement, in the thermodynamic limit N — > oo, defined 
as 



Sa 
N ' 



(4) 



We will show below that a nonzero entanglement exists 
in the thermodynamic limit, which will give us a han- 
dle for examining the entanglement of quantum phase 
transitions. In contrast, the bipartite entanglement be- 
tween two spacial blocks does not have the thermody- 
namic limit, as the block entropy does not scale with 
the system size 0, Q. A quantum phase transition is 
usually accompanied by a singular behavior of a thermo- 
dynamic observable as the coupling strength is tuned to 
a critical value, similar to the thermal phase transition 
where a thermodynamic potential exhibits a singularity 
as the temperature is tuned to a critical value. We expect 
the multi-species entanglement described above is capa- 
ble of highlighting the quantum transition in the ther- 
modynamic limit. Below, we will study the macroscopic 
entanglement in the context of the Ising model in the 
presence of a transverse magnetic field, which is a proto- 
type exactly-solved model of quantum phase transitions. 

Let us consider N lattice sites occupied by a spin-1/2 
species, whose interactions arc described by the Hamil- 
tonian given by, 



% = — J Y^cr; rr 



- h 



(5) 



Here, of, of are pauli operators acting on the local 
hilbert space at site i. Let | f), | I) be the of eigenstates 
with of = ±1 respectively. The system exhibits long- 
ranged order, either ferromagnetic (for J > 0) or antifcr- 
romagnetic (for J < 0) when h — 0. The ground state is a 
direct product of of eigenstates, |±) = | f)±| L)/v2, the 
ferromagnetic state having all sites in |+) state where as 
the antiferromagnetic state has alternating sites in |+) 
and |— ) states. As the transverse field is switched on, 
quantum fluctuations induce excitations, and the ground 
state is a coherent superposition of h — ground and ex- 
cited states, with the weights depending on the coupling 
strength. The situation is similar to inducing transitions 
to excited states, through coupling to a thermal bath, 
except that the state of the system in that case would be 
an incoherent mixture of states, with the temperature de- 
termining the weights. For \J\ > h, the system exhibits 
long-ranged order, and for | J| < h there is no long-ranged 
order. The system exhibits a quantum critical behavior 
at | J | = h, with a vanishing excitation energy gap@ . 
The concurrence measure of entanglement between near- 
est neighbor sites is shown to be singular as a function 
of the coupling strength, the first derivative showing a 
logarthmic divergence^. We will show below that near 
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FIG. 1: The integrated density of eigenvalues as a function 
of the eigenvalue location, for a few values of the interaction 
strength x = J/h. 



the quantum critical points x = J/h = ±1, the entangle- 
ment e(x) between up and down spins exhibits singular 
behavior, with its derivative showing a jump discontinu- 
ity and the its second derivative diverges as ||jc| — 1| . 

The above Hamiltonian has been exactly solved for all 
eigenstates by mapping the spins to fermions through 
the Jordan- Wigner transformation [|| Q. followed by a 
Fourier transformation to momentum basis, followed by a 
Bogoliubov or quadratic-form diagonalization. The first 
step is to go from the spin variables of , of = + of 
to fermion variables rii , cj , a the fermion number opera- 
tor, the creation and annihilation operators, through the 
Jordan- Wigner transform given as, 



of = 2ni — 1, a/ 



e i7r £i=i™3 C J 



(6) 



The presence of the phase-factor operators in the above 
makes the multi-spin correlation functions very difficult 
to calculate, however for computing the macroscopic en- 
tanglement they do not pose a great difficulty. In terms 
of the fermion states, the presence of a fermion at a site 
(i.e. rii = 1) implies that site is occupied by an up spin, 
and an absence of a fermion implies the site is occupied 
by a down spin. Now, let us define the momentum-basis 
fermion operator, given by 



—y 



c\ e 



-iql 



The allowed values of q for the case of periodic boundary 
conditions and N even, are q 



I 7T I 37T 



± 



(JV-I)tt 



for 



the sectors with even number of fermions Np. Similarly, 
for Np odd, the allowed values are q — 0,±t5,±4£... ± 

n ' 7r ' ^ ne Hamiltonian becomes uncoupled in terms 
of different \q\ values, and for each q > 0, a simple di- 
agonalization from the basis |0), \(f> q ) = cJcljO), | — q) = 
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FIG. 2: The entanglement, and its first derivative, between the up and down spin particles is plotted as a function of the 
coupling strength, for N = 10, 20, and in the thermodynamic limit N — ► oo. 



c! g |0), \ q) = cj|0) yields the eigenvalues and eigenstates 
of the Hamiltonian. Of particular interest is the ground 
state, is given as a direct product 0], 



iG)=nKio)+ 



g>0 



The amplitudes depend on the coupling strength, and 
thus determine the entanglement in the state. We have, 



h + J cos q 



J 2 + 2 Jh cos q 



)• 



(9) 



and the other amplitude is given by \b q \ 2 = 1 — |a 9 | 2 . 

The above state is a superposition of many-particle 
momentum basis states, the amplitude for a given mo- 
mentum basis state would be a product of wave func- 
tion amplitudes, either a q or b q appearing for each value 
of q, corresponding to either the q state is occupied by 
fcrmions (up spins) or the unoccupied by fermions (or 
equivalently occupied by holes, i.e. down spins). Thus 
each component of the above state can be labeled by u (as 
in Eq. 2), a set of q values of the occupied states of the 
fermion (up-spin particles) , owing to the exclusion prop- 
erty, the complement of the set u would be the occupied 
q values for the unoccupied (down-spin particles). Thus, 
after tracing over the down-spin degrees of freedom, the 
up-spin state is obtained as, 

pt=n a ?i°)(°i+E^ n 4im--^-.o)<o,o..<v-oi+.. 

q q'^q 

(10) 

Here, the first component is the contribution of the vac- 
uum state, and the second series term is the contribution 
from the two-particle sector and so on. The structure of 
the reduced density matrix is similar even for the maxi- 
mum energy state (each of the amplitudes a q is replaced 
by the amplitude b q and vice versa), and hence the same 
quantum critical behavior from the view point of the two- 
species macroscopic entanglement. 



The von Neumann entropy can be shown to be a sum 
of independent contributions of q modes, the individual 
mode contribution being equal to the Shannon binary 
entropy H(pi) = -p t logp 4 - (1 - pi) log(l - pi), where 



(8) the eigenvalue is given by p. 



As the coupling 



strength x — J/h is varied, the band of eigenvalues Pi 
and the band of eigenvalues 1-pi spread in width. For 
|x| < 1, as qi is varied, \a q \ 2 is bounded well below the 
value of 1/2, and similarly \b q \ 2 well above 1/2, with a 
gap A (a;) between the two bands of eigenvalues. The 
gap decreases as x is increased, and finally the gap van- 
ishing at | a; | = 1. For |at| > 1, the values of |a 9 | 2 (thus 
\b q \ 2 ) span the full range. A gap in the individual mode 
eigenvalue would imply a gap between the largest and the 
second largest eigenvalues of the reduced density matrix 
shown above. From the functional relation of the wave 
functions given in Eq. 10 we can work out the gap, and 
we have, 



A(x) = \fl - x 2 0(1 - \x\), 



(11) 



which gives a gap exponent of 1/2, i.e. A ~ (1 — |a; | ) x / 2 
as \x\ 1. In contrast, the particle- hole excitation 
gap from the ground state to the excited state vanishes 
linear lyQ. 

The macroscopic entanglement (from Eq.5) between 
up and down spin particles for a given coupling strength 
x = J/h, is given by 



e{x) 



i N/2 r 1 



dp g(p, x) lo{ 



I-P 

P 



(12) 



We have converted the sum into an integral by introduc- 
ing the integrated density of eigenvalues, i.e. the number 
of eigenvalues pi below a given value of p, defined by 



(13) 
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Here, Pi(x) = (1 — Q)/2 is the eigenvalue location, Q = 
\l+x cos q\/\/l + x 2 + 2a;cosg. Thus, we have transferee! 
the dependence on the coupling strength and the system 
size to the density of eigenvalues. The integrated density 
of eigenvalues can be calculated in the thermodynamic 
limit, the details will be published elsewhere. Using p = 
(1-0/2, we have g( P ,x) = \ + 6(()9(( 2 + x 2 - 1)(*_ - 
<f> + )/27r, where 



C0S z ± = z±±f±mEMVe±ZEJi. (14) 

X 

The integrated density of eigenvalues is plotted in Fig. 1 
for a few values of the coupling strength x. For p > 1/2, 
it saturates to the value of 1/2 for all values of x, and 
for |a;| > 1, the function g(p,x) has no dependence on x\ 
The entanglement now is given by, 



dClogi±|(^ 



-$+)0(x 2 +C 2 -l). (15) 



We have plotted in Fig. 2 the entanglement and its 
first derivative as a function of x = J/h from the sum 
of Eq.12 for N = 10, 20, along with that of the infinite 
chain from equation 15. The entanglement reaches its 
maximum value of about s w 0.19 near |x| = 1. For 
the finite-size cases, the entanglement increases a little 
after |a;| = 1, and quickly falls to this value. Contrast 
this, with the entanglement of log 2 for the ground state 
with h = 0, implying that the limit of h — > is not the 
same as h = 0. Even for a finite size cluster of N = 10 
sites, one can see in the figure a clear signature of the 
jump discontinuity near the quantum critical points. The 
value of the macroscopic entanglement as x — > oo can be 
calculated, as $± — > ±C, we get e(x — > oo) w 0.19. 



The first derivative shows a discontinuity at |a;| = 1, 
which corresponding to the quantum critical point, and 
the second derivative shows a divergence. Since the be- 
havior of the entanglement depends on the functions $±, 
and de/dx depends on $' ± , we can expand these functions 
near the quantum critical points. We have the behavior 
near the critical point x =!, $' + w —y/l — ( 2 /( which 
does not exhibit any discontinuity, and exhibits a 
discontinuity at x = 1, we have 



dx 



1-X y/W^ 

|l-a;| C 



(16) 



Thus, for x — > 1+, s' = 0, but the derivative is nonzero 
as x — > 1~ . The discontinuity 5s' of the derivative of the 
entanglement at the quantum critical point is given by 



6s' (x = 1) w 0.28. 



(17) 



Now, the divergence of d 2 e/dx 2 depends on d 2 & ± /dx 2 . 
Near the quantum critical point, $ + has a finite second 
derivative, but <£>_ shows a singularity, as shown by, 



Q 2 <i>_ 



l v/w 5 



dx 2 ~\x-i\ ( ■ (18) 

Hence, the second derivative of the entanglement diverges 
near x = 1, as s" <~ \x — 1| _1 - 

In conclusion, we have shown that the multi-species 
macroscopic entanglement, nonzero in the thermody- 
namic limit, exhibits singular behaviour near quantum 
phase transitions. It will be interesting to see if the dy- 
namics starting from a given initial state can also cap- 
ture signature of critical behaviour, and whether finite- 
temperature behavior can be explored. 
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